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' We compare different methods of obtaining the neutrinoless double beta decay nuclear 

matrix elements (NME). On the example of '^^Ge we use the NME to calculate the 
Majorana neutrino transition magnetic moments, generated through particle— sparticle 
^ , ij-parity violating loop diagrams whithin the minimal supersymmetric standard model. 

1. Introduction 

The neutrinoless double beta (Oj^2/3) decay is a hypothetical nuclear process, in 
' which two simultaneous beta decays occur in one nucleus, whereas some exotic 

04 . mechanism prevents the neutrinos to be emitted. It would result of course in the 

lepton number violation (AL = 2) and therefore it is forbidden by the standard 
theory. However, in certain exotic models of physics beyond the standard model, 
such processes are allowed and sought in a number of experiments. The observation 



X 



$H ' of a Oz/2/3 decay will open a completely new area of physics, which by far is only 

^ a speculation. 

The theory of this exotic decay involves two parts: firstly, as this is a nuclear 
process it requires a careful calculation of the nuclear matrix elements; secondly, 
an exotic mechanism must be proposed, which would suppress the emission of the 
neutrinos. The Gv2fi half-life may therefore be expressed in a factorized form as 

(T""-)-! -G°nM°'^n(mee)|', (1) 

where G"'^ is the exactly calculable phase-space factor, M*^" is the NME, and (rriee) 
is the so-called effective neutrino mass, which represents the non-standard interac- 
tions involved. We are not going to discuss the latter part here. 

In this communication, we are going to compare different methods of obtaining 
the nuclear matrix elements (NME) for the 0i^2/3 decay. As an application, we will 
use them to calculate the Majorana neutrino transition magnetic moments within 
the minimal supersymmetric standard model with broken i?-parity. 
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2. Calculations of the NME: random phase approximation 

The random phase approximation (RPA) methods in their simplest form base on 
the assumption, that the ground state is expressed as a BCS vacuum. Then, a tran- 
sition amplitude between the 0"^ BCS vacuum of an even-even nucleus and excited 
states of the neighboring odd-odd nucleus is calculated, representing the transition 
between the mother nucleus and an intermediate one. A similar amplitude repre- 
sents the transition between the intermediate and the daughter nucleus. These two 
parts are then summed over possible excited states of the intermediate virtual nu- 
cleus, which are expressed as simple harmonic oscillations above its BCS ground 
state. 

Such simplified picture does not give satisfying results and many variation of the 
original method appear. These involve the use of quasiparticles (QRPA, pnQRPA), 
renormalized QRPA (RQRPA), selfconsistent RQRPA (SRQRPA) and others. In 
QRPA a quasiboson approximation is used, in which the commutation relations 
for fermions are replaced by bosonic ones, and this procedure seems to be working 
well for small harmonic excitations, but breaks down when one wants to take into 
account higher excitations as well. In RQRPA the commutators are approximated 
in a different way, which introduces the new factor as a normalization constant. 

What is more, different schemes of short-range correlations are used, among 
which the Jastrow and the unitary correlation operator method (UCOM) are most 
common. Also one may adjust other parameters, such as the nuclear potential 
(Woods-Saxon, effective mean field etc.), and the strength of the particle-particle 
(gpp) and particle-hole (gph) interactions. 

All these give a quite big diversity in the outcome of different calculations. At 
present one may summarize the results of the Jyvaskyla and the Tiibingen groups,^ 
for ™Ge as 3.33 < AfO" < 6.64. 



3. Calculations of the NME: nuclear shell model 

The nuclear shell model (NSM) and its recently announced large-scale version 
(LSSM) is meant in principle to give exact results. The approach is straightfor- 
ward and involves defining a valence space, deriving an effective interaction from 
the Hamiltonian, building, and then finally diagonalizing the Hamiltonian matrix. 
The first and most obvious obstacle here is, however, the dimension of the matrix 
in question, which is proportional to 




where (d^) is the dimension of the proton (neutron) subshell, and p (n) is the 
number of valence protons (neutrons). This feature makes it practically impossible 
to include in the valence space all important single particle orbits; eg. the spin- 
orbit partners are usualy neglected. Also in the actual calculations the many-body 
problem is reduced to two-body problem within some mean-field approach. 
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Table 1. Neutrinoless double beta decay nuclear matrix element 
M'^" for ''^Ge calculated using three different approaches. The 
last column shows the corresponding upper limit on the effective 
neutrino mass (see text for details). 



method 




ranges 


central value 


Km 


ee)\ < 


(R)QRPA 

LSSM 

IBM 


3.33 
2.22 
4.64 


- 6.64 

- 2.81 

- 5.46 


4.985 ± 1.655 
2.515 ± 0.295 
5.050 ± 0.414 


0.22 - 
0.51 - 
0.26 - 


0.43 eV 
0.65 eV 
0.31 eV 



The most advanced large scale shell model calculations for the Ge nucleus 
yield at present M"'^ = 2.22 (without higher order contributions to the nuclear 
current this value increases to 2.58). A little bit older calculations which used the 
UCOM short-range correlations gave M°'' = 2.81.^ 



4. Calculations of the NME: interacting boson model 

The interacting boson model (IBM) has not been used for this type of calculations 
before, but it turned out to give results which surprisingly well agree with the 
(R)QRPA calculations. In this method, a fermionic two-body matrix element is 
used to obtain the fermions transition operator in second quantization. The matrix 
elements of this operator are then evaluated in the general seniority scheme. At 
the end, the fermionic operator is mapped into bosonic transition operator, whose 
matrix elements are evaluated using realistic IBM wave function. 

This method, introduced by Barea and lachello in Ref.^, yields for the ^^Ge 
nucleus 4.64 < M"'' < 5.46, which is in very good agreement with the central 
values obtained within the (R)QRPA method. 

We summarize the various results for the ^^Ge nucleus in TablTJ In the last 
column we have calculated the effective neutrino mass using the result of the 
Heidelberg-Moscow experiment, T^J^(J^Ge) < 1.9 x 10^^ y,-'^ and a common phase- 
space factor G^" = 2.55 x 10^26 y cV^. 



5. Majorana neutrino transition magnetic moments in i2-parity 
violating supersymmetry 

The Majorana magnetic moment acts between and v^j^ chiral components of 
Majorana neutrinos, assuming a standard gauge theory with only left-handed neu- 
trinos. As a consequence, it violates the total lepton number by two units AL = 2. 
The effective Hamiltonian iJoff of this interaction takes the form 

^^cff = ^fJ-zji^zLCr^'^'^jLFlSa + h-C, (3) 

where F is the electromagnetic field strength tensor, and a"^ is defined by means 
of the Pauli matrices, a"^ — ^[c",cr^ . We work in the minimal supersymmetric 
standard model with broken i?-parity^ In this setting it is possible to consider 
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Fig. 1. Processes leading to Ma- 
jorana neutrino transition mag- 
netic moments. 



processes, in which the Majorana neutrino, despite being electrically neutral, may 
effectively interact with an external photon. The amplitude of such processes is 
governed by the magnetic moments of the neutrino. We consider interaction between 
two Majorana neutrinos, in which the effective interaction vertex is expanded into 
a particle-sparticle loop, as depicted on FiglT] 

The quark-squark loop contribution, including possible d-quark mixing {V is 
the CKM matrix), is given by the following formula:^ 

= (1 - E E ^^^^^^^^ - ^ j)| (4) 

jkl (a ) 

where w is the loop integral, which arises from the integration over virtual momenta 
of particles inside the loop, 

„ sin 26**^ / x-!^ Inx-'J' — xi^ + \ , ^^ 

and we have denoted the ratio of the quark and squark masses squared by x-'^ = 
m?j^j/m?-^. Here 9^ is the fc-th squark mixing angle. The contribution from the 
slepton-lepton loop does not contain summation over three quark colors, and we 
do not include the very weak mixing in the charged lepton sector. The formula 
reads thereforJ^l 

= (1 - '^-o^^ E [:^r^.) 

where the loop integral is identical as previously with the exception that x''^ — > 
yl = m'jj/m'j^, and 6'^ 4>^ , the slepton mixing angle. 

In order to get numerical results, we constrain the model by imposing grand 
unification at very high energy scale, Mqut ~ 1.2 x 10^^ GeV, thus starting from few 
free parameters only. These are the common mass of scalars mo, the common mass of 
fermions fTii/2, common Yukawa coupling constants unification factor Aq, the ratio 
of the Higgs vacuum expectation values tan /3, and the sign of the bilinear up- and 
down-type Higgs coupling sgn(/i). Next, the renormalization group equations are 
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LSSM: 0.816 - 






LSSM: 0.789 | 


(R)QRPA: 0.668 


1 IBM: 0.672 






(R)QRPA: 0.571 






^ IBM: 0.555 









Fig. 2. Ranges for the maximal values of the off- 
diagonal elements in the neutrino mass matrix, 
coming from different 0i'2[5 nuclear matrix ele- 
ments. 



Table 2. Majorana neutrino transition magnetic moments /i^^^ in /ib for 
GUT parameters: Aq = 100 GeV, mo = mj/a = 150 GeV, tan/3 = 19, 
fi > 0. Ranges correspond to the spread in NME obtained using different 
methods. 

ij LSSM (R)QRPA IBM 

lepton-slepton loop mechanism 
e^i,eT (1.33,1.33) 10-15 (8.46, 13.3) 10-i« (9.99,11.9)10-16 
fiT (1.23,1.35) 10-15 (7.02,11.6) lO-i'' (8.28,9.86) lO-^*^ 

quark-squark loop mechanism (without d-quarks mixing) 
efi,eT (9.49,9.49) 10-^^ (6.04,9.49) lO-i'^ (7.13,8.50) lO-i'^ 
IMT (8.69,9.52) lO-i'^ (4.95,8.22) lO-i'^ (5.84,6.96) lO-i'^ 

quark-squark loop mechanism (with d-quarks mixing) 
e/i,eT (8.22,8.22) lO-l'^ (5.24,8.22) lO-l'^ (6.18,7.36) lO"!'^ 
IMT (7.24,7.93) lO-i'^ (4.13,6.85) IQ-'^'^ (4.87,5.80) lO"" 



used to derive the values of the couphngs and mass parameters at low energy scales. 
The i?-parity violating trilinear couplings A and A', which cannot be derived from 
the GUT constraints, are assessed from the experimental neutrino mass matrix. 
This procedure allows to calculate the transition magnetic moments (see for 
details). 

FigO presents how the off-diagonal elements in the neutrino mass matrix de- 
pend on the nuclear matrix element used. Each bar shows the range for a given 
method, with the number being its central value. One sees immediately, that the 
(R)QRPA and IBM methods agree very well with each other, while the shell model 
calculations give substantially higher results. This trend is of course preserved also 
in the magnetic moment data, which are shown in Tabsl2]and|3l We have included 
calculations for two cases, one with 'small' GUT scale unification parameters, the 
other with 'large' ones. For each case three possibilities were considered, namely 
the lepton-slcpton loop, and the quark-squark loop with or without additional d- 
quark mixing. One sees that the differences in the magnetic moments for various 
nuclear matrix elements are negligible, and that the ranges often overlap with each 
other. As expected, the LSSM gives the highest results, but the difference between 
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Table 3. Majorana neutrino transition magnetic moments //i/-^ in fiB for 
GUT parameters: Aq = 500 GeV, mo = mi/2 = 1000 GeV, tan /3 = 19, 
/I > 0. Ranges correspond to the spread in NME obtained using different 

methods. 



LSSM 



(R)QRPA 



IBM 



Icpton-slepton loop mechanism 



(3.48,3.48) 10"" (2.22,3.48) 10"" (2.62,3.12)10"" 
(3.26,3.57) 10"" (1.86,3.08) 10"" (2.19,2.61)10"" 



HT 



quark-squark loop mechanism (without d-quarks mixing) 

(2.53.2.53) 10^1** (1.61,2.53)10^1** (1.90, 2.27) 10"!^ 

(2.32.2.54) 10"!* (1.32,2.19)10^1** (1.56,1.86) 10"!* 



e/i, er 



quark-squark loop mechanism (with d-quarks mixing) 
(2.15,2.15)10^1** (1.37,2.15)10^1** (1.62,1.93) 10^1** 
(2.03,2.22) 10"!* (1.16,1.92)10"!* (1.36,1.62) 10"!* 



for example the upper (R)QRPA value and the lower LSSM one is roughly of the 
order of 10%. Of course, the predictions on the half-life of the Oz/2^ decay, or the 
effective neutrino mass may differ in a much more serious manner. In such a case 
the (R)QRPA and IBM methods are worth recommendation. 
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